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Abstract 

In the framework of the minimal supersymmetric standard model and the 
i?-parity breaking model, we investigate various production processes of the 
supersymmetric partner at HERA energies. Our emphasis is paid upon the 
scalar top quark, the partner of top quark, characterized by its lighter mass 
than the top quark and other scalar quarks in a model. We propose experi- 
mentally feasible approaches to search for clean signals of the stop from either 
its production or decay processes. 
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1 Introduction 



It is very amazing that the Standard Model (SM) has been so successful [|l|]. Discoveries of 
and Z° and high level precision electroweak tests by LEP have been vindicating the 
SM over and over again. The discovery of the top quark at TEVATRON |Q] seems to be 
its culmination. The theoretical prediction agrees very well with experiments. However, 
the problems of and R c may shed a light on the fate of the SM ||. As far as the 
first approximation is concerned there is no need to go beyond the SM. However, the SM 
cannot be a fundamental theory in any sense. Actually, there are a large number of free 
parameters, the arbitrariness of particle masses and mixing angles, and the lack of any 
explanation for the replication of generations and so on. Moreover, it is known that the 
gauge hierarchy problem exists in the SM. Theorists have so far attempted to find a way to 
go beyond the SM via various standpoints : unification (|], technicolor ||, supersymmetry 
(SUSY) § or superstring §. 

Among them the SUSY seems to be the most reasonable candidate to our expectation. 
A number of theoretical and phenomenological reasons make low energy SUSY attractive 
with respect to its alternative. The supersymmetry is a symmetry between fermions and 
bosons. In the SUSY models, the quadratic divergence can be cancelled out owing to both 
contributions of boson and fermion loops. This is a key point for the fact that the hierarchy 
problem can be solved in the SUSY models. However, we should note that fermions and 
bosons, in other words, the ordinary particles and their SUSY partners (sparticles), must 
be degenerate in mass in the exact SUSY limit. Obviously there appears to be no evidence 
in nature for such a situation. Therefore, in order to apply SUSY to particle physics, we 
must consider models in which SUSY is broken. In this case masses of all sparticles must 
be less than about 1 TeV in order to solve the hierarchy problem. At first sight, it seems 
to be unnatural that we need a large number of new particles, sparticles, undiscovered yet. 
However, such new particles play an essential role in some phenomenologically favourable 
properties in the model. The most impressive evidence in favour of SUSY may be the 
unification of gauge couplings in SUSY Grand Unified Theories (GUTs) ||. The dark 
matter issue || in astrophysics also encourages SUSY proponents. Most of us consider 
the "SUSY world" as a plausible scenario for future particle physics. 

HERA at DESY is the world first electron-proton collider whose ep center of mass en- 
ergy is 314 GeV produced by 30 GeV electron (positron) and 820 GeV proton. The design 
luminosity is 1.5 xlO 31 cm~ 2 s . Polarized electron and positron beams are available. Two 
experimental groups HI and ZEUS are now engaged in various experiments. Needless to 
say, HERA is expected to be the powerful machine to serve us utmost information on the 
nucleon structure. In addition, HERA has a large potentiality to discover signatures of 
the new physics beyond the SM. This is because of higher energies available than LEP and 
cleaner events than hadron colliders such as TEVATRON. Also the unique signature at 
the ep collider which is not obtained from electron-positron colliders nor hadron colliders 
reveals us a striking feature of new physics. 

The purpose of the present paper is to review the search for sparticles at the ep 
collider HERA. Our emphasis will be paid upon the scalar top quark (stop), the partner 
of top quark, characterized by its lighter mass than the top quark and other scalar quarks 
(squarks) in a model described later [[h], |llj. The existence of the light stop could give 
a clue to the issue on [12]. The expected mass of the stop is clearly within the reach 
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of HERA. Particularly, the squarks could be singly produced in ep collisions |13|, 14, [Tf| 
in the i?-parity breaking model (RBM) (l(| because of the existence of the electron- 
quark-squark couplings. Clearly HERA is the best machine to search for the stop in 
such models, because we expect the remarkable peak structure in the Bjorken parameter 
x distribution. We calculate cross sections together with detectable signals expected to 
observe at HERA. Although we have studied production processes for selectron, sneutrino 



or photino at HERA [17] we do not enter their details because of their discovery potential 
to be not so large. 

The organization of the paper is as follows. In Sec. 2 we develop the theoretical 
framework of Minimal SUSY Standard Model (MSSM) and the origins of the light stop 
mass is discussed. Decay modes of the stop with mass less than other squarks and gluino 
are also exploited. Experimental status of SUSY particle search is extensively presented 
in Sec. 3. Not only status on HERA, but also data on the search from LEP, TEVATRON 
and TRISTAN are concisely summarized. Various analyses on the sparticle production 
and their expected signatures from decay processes at HERA are reviewed in Sec. 4. 
Stop production by boson-gluon fusion and i?-parity breaking single stop production are 
presented in detail. The concluding remarks are given in the final section 5. 



2 Theoretical framework of MSSM 
2.1 Particle content 

The Minimal SUSY Standard Model (MSSM) includes the minimal particle content. That 
is, there should be a new particle (sparticle) for each known particle in the SM and the 
two Higgs doublets. The additional Higgs doublet must be included in order to produce 
masses of up- and down- type quarks and to realize the chiral anomaly cancellation ||. 
The list of particles in the model is shown in Table I (weak eigenstates) . 

For quarks and charged leptons there should exist two scalar partners per species. 
These scalars are called squarks and sleptons, generically sfermions. Before SU(2)xU(l) 
breaking, the left- and right-sfermions do not mix since they have different SU(2)xU(l) 
quantum numbers. After the breaking, however, they can mix each other. Actually, this 
mixing effect is substantial only for the third-generation sfermions, especially for the stops 
(superpartners of the top quark). The left-right mixing of the stops will be discussed in 
the following subsection. 

The new fermions are either the superpartner of spin-1 gauge bosons (gauginos) or 
that of spin-0 bosons (higgsinos). They mix each other when SU(2)xU(l) symmetry is 
broken. The mass eigenstates (inos) are usually some complicated mixture of gauginos and 
higgsinos. Electrically neutral-inos and charged-inos are called neutralinos (k = 1 ~ 4) 
and charginos Wi (i = 1,2), respectively. Gluinos g are free from the mixing since the 
color SU(3) is not broken. 

We need two Higgs doublets in the MSSM [18|. We have five physical Higgs bosons 
which remain after SU(2)xU(l) breaking. They are two CP even (h°, H°) and one CP 
odd (^4°) neutral scalars and remaining one charged scalar (H^). 
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2.2 Basic parameters of MSSM 



The standard model has 18 fundamental parameters to be determined by experiments. 
The MSSM has a somewhat larger number of parameters. They are classified as (i) gauge 
couplings, (ii) superpotential parameters and (iii) soft-breaking parameters. 

2.2.1 gauge couplings 

The three gauge coupling parameters corresponding to SU(3), SU(2) and U(l) gauge 
groups are the same as in the SM. All gauge interactions are governed by these cou- 
plings. They determine the fermion - sfermion - gaugino interactions and four-point scalar 
interactions as well as ordinary fermion - fermion - gauge-boson interactions. 

2.2.2 superpotential parameters and R-parity 

In the MSSM ordinary Yukawa interactions are generalized to the terms in superpoten- 
tial W((f>), where <fi denote arbitrary chiral superfields. Renormalizability restricts the 
functional form of the superpotential to 

W(4>) = mij<j>i<j)j + \jk<i>i<i>j<i>k (!) 

The parameters rriij and A^-fe are further constrained by the gauge symmetry and some 
discrete symmetries. 

A well-known discrete, multiplicative symmetry is the .R-parity defined by 

R=i _f(B-L) + 2S^ (2) 

where B, L and S denote the baryon number, the lepton number and the spin, respec- 
tively. This formula implies that all ordinary SM particles have even .R-parity, whereas 
the corresponding superpartners have odd .R-parity. Usually we impose the B — L conser- 
vation on the MSSM and then the MSSM possesses the imparity invariance. In this case 
the superpotential has the form ; 

W = nH x H 2 + hBi&B + f d HiD°Q + f u H 2 U c Q, (3) 

where we have used the superfield notation in Table I. Here parameters /'s correspond to 
usual Yukawa couplings and fi is a new supersymmetric parameter which has no counter- 
part in the SM. 

The .R-parity conservation in scattering and decay processes has a crucial impact on 
the SUSY phenomenology at high energy colliders. First, the SUSY particles must be 
produced in pairs. Second, the lightest SUSY particle (LSP) is absolutely stable. On 
the other hand, the cosmological constraints tell us that the LSP should be electrically 
and also color neutral. Consequently, it is weakly interacting in ordinary matter. The 
heavy unstable SUSY particles must finally decay into the LSP. Therefore the canonical 
signature for the .R-parity conserving SUSY models at collider experiments is the large 
missing (transverse) energy. 
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It should be emphasized that the imparity conservation is not automatically imposed in 
the MSSM. In fact the i?-parity breaking superpotential is allowed by the supersymmetry 
as well as the gauge symmetry ; 

Vfy= \ ijk Ul 3 E% + X' ijk LiQjD c k + Xl jk U c iD^D c k . (4) 

The first two terms violate L and the last term violates B. If we want to explain some 
unresolved problems such as (i) the cosmic baryon number violation, (ii) the origin of the 
masses and the magnetic moments of neutrinos and (iii) some interesting rare processes 
in terms of the L and/or B violation, the -R-parity breaking terms must be incorporated 
in the MSSM. 

In such a model the signatures for SUSY particles observed at collider experiments 
should be very different from the caconical one. We will discuss the typical/? signals at 
HERA in Sec.4. 



2.2.3 soft-breaking parameters 

SUSY is not an exact symmetry of nature since our world is not manifestly supersymmetric. 
In the MSSM the SUSY breaking is induced by the soft-SUSY breaking terms, which do not 
introduce quadratic divergences. Hence the solution of the naturalness problem remains 
intact. There are four types of soft breaking terms ; (i) gaugino masses Mj (i = 1 ~ 3), 
(ii) masses for the sfermions m/, (iii) trilinear term Af and (iv) three scalar Higgs mass 
terms. These three mass parameters can be re-expressed in terms of two Higgs vacuum 
expectation values, v\ and v%, and one physical Higgs mass. Here v\ and V2 respectively 
denote the vacuum expectation values of the Higgs field coupled to d-type and n-type 
quarks. v\ + v\ = (246 GeV) 2 is determined from the experimentally measured VF-boson 
mass, while the ratio 

tan/3 = — (5) 

Vl 

is a free parameter of the model. 

A generally accepted assumption is that all the three gaugino mass parameters Mj 
are equal at some grand unification scale Mx- Then the gaugino mass parameters can 
be expressed in terms of one of them, for instance, The other two gaugino mass 

parameters are given by 

Mi = ^ tan 2 9 W M 2 (6) 

M 3 = -sin 2 8 W M 2 , (7) 
a 

where a and aj, denote the QED and the QCD coupling constants, respectively. 



2.3 Scalar top in MSSM 
2.3.1 left-right mixing of stops 

In the framework of the MSSM Q, the stop mass matrix in the (t^,, tn) basis is expressed 
by 

/ mi- a+mt \ 

M 7=[ tL V • (8) 
y a t m t J 



5 



where mt is the top mass. The SUSY mass parameters and at are parametrized in 



the following way [19] : 



to 2 



to 2 - 
at 



m 



q 3 + % cos 2(3 
2 

c/3 + 3 
A t + (i cot /3 



sin 2 ) + m t- 



mfj + -m 2 z cos 2(3 sin 2 + m 2 . 



(9) 

(10) 
(11) 



The soft breaking masses of the third generation doublet mg 3 and the up-type singlet mu 3 
squarks are related to those of the first (and second) generation squarks as 



m Q 3 



nr 



21, 



(12) 
(13) 



where / is a function proportional to the top quark Yukawa coupling at and is determined 



by the renormalization group equations (RGEs) [2C] in the minimal supergravity GUT 



(MSGUT). Throughout this paper we adopt the notation of Ref.|21]. 



There are two origins for lightness of the stop compared to the other squarks and 
sleptons, i) smallness of the diagonal soft masses to 2 - and to 2 - and ii) the left-right stop 

*L tR 

mixing. Both effects are originated from the large Yukawa interaction of the top. The 
origin i) can easily be seen from Eqs.@~(|l3|). The diagonal mass parameters to 2 - and 

to 2 - in Eq.dq) have possibly small values owing to the negative large contributions of I 

tR ^ 

proportional to at in Eqs.(p^) and (|r3|). It should be noted that this contribution is also 
important in the radiative SU(2)xU(l) breaking in the MSGUT. The Higgs mass squared 



has an expression similar to Eqs.(|r2|) and (13) ; 



TO 



H 2 



3/, 



(14) 



where to 2 denotes the soft breaking mass of the first generation doublet slepton. The 
large contribution of I enables to 2 ^ to become negative at an appropriate weak energy 
scale. In order to see another origin ii) we should diagonalize the mass matrix (||). The 
mass eigenvalues are obtained by 



TO- 



2 2 ( ( 

+ to- =f ( ( 



t L 



tR 



t L 



mi/ + {2atmtf" ' " 



(15) 



and the corresponding mass eigenstates are expressed by 

fh\ _ (t L cos 9 t - t R sin 9 t 
\t2 J \Jl sin6> t + t R cos9 t/ 

where 9t denotes the mixing angle of stops : 



(16) 



tan#/ 



a t m t 



2 

tR 



2 



(17) 



We see that if the SUSY mass parameters and the top mass are of the same order of 
magnitude, small is possible owing to the cancellation in the expression (15) [10, 11]. 
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After the mass diagonalization we can obtain the interaction Lagrangian for the mass 
eigenstate t\. We note, in particular, that the stop coupling to the Z-boson (t\t\Z) 
depends sensitively on the mixing angle Ot- More specifically, it is proportional to 



= | sin 2 0w - i cos 2 t . UN 



Note that for the special value of #t~0.98, the Z-boson coupling completely vanishes [22] 



2.3.2 R-parity breaking interactions of stop 

Interesting properties of the lighter stop ti described in the last sub-section are not mod- 
ified even in the i?-parity breaking models. In this case the stop could have interactions 
with the ordinary leptons and/or quraks via the i?-breaking Lagrangian ; 

L = X' i3j cos OthdjPLli + h.c, (19) 

which is originated from the superpotential (Q). In particular the interaction with the 
electron ; 

L = A' 131 cos 6 t {tidP L e + t\eP R d) (20) 
will be most suitable for the ep collider experiments at HERA because the stop will be 



produced in the s-channel in e-q sub-processes [13]. Note that the stop can not couple to 
any neutrinos via /^-breaking interactions. This is a unique property of the stop which 
could be useful for us to distinguish the stop from some leptoquarks. 



2.3.3 decay modes of stop 

Here we examine the decay modes of the stop. In the MSSM, the stop lighter than the 
other squarks and gluino can decay into the various final states : 



tz k 


(a) 


bW t 


(b) 


blv 


(c) 


bvl 


(d) 


bWZ k 


(e) 


bfJZk 


(f) 


cZ%, 


(g) 


ed, 


00 



where Zk (k = 1 ~ 4), Wi(i = 1,2), v and I, respectively, denote the neutralino, the 
chargino, the sneutrino and the charged slepton. (a) ~ (g) are the i?-parity conserving 
decay modes, while (h) is realized by the RB couplings (|20|). 

If we consider the stop with mass small enough in the R conserving case, the first five 
decay modes (a) to (e) are kinematically forbidden due to the observed top mass mt~175 
GeV p as well as the model independent lower mass bounds for sparticles ; m~ ~ 45 
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GeV, m~ ~ 45 GeV and m~ ~ 40 GeV. So (f) and (g) survive. Hikasa and Kobayashi 



[11] have shown that the one-loop mode t\ — ► cZ\ (g) dominates over the four-body mode 
t\ — > bf f Z\ (f). So we can conclude that such a light stop will decay into the charm 
quark jet plus the missing momentum taken away by the neutralino with almost 100% 
branching ratio. On the other hand, if we consider the RB coupling A' 131 > 0.01, which 
roughly corresponds to the coupling strength detectable at HERA, the decay modes (c) 
to (g) are negligible due to their large power of a arising from multiparticle final state or 
one loop contribution. So there the two body modes (a), (b) and (h) survive in this case. 



3 Experimantal status of SUSY particle search 

Most of the theoretical models assume -R-parity conservation, which has important physical 
implication: (1) SUSY particles must be produced in pairs, (2) heavy SUSY particles decay 
to lighter SUSY particles, and (3) the LSP is stable. The phenomenology of imparity 



violating interactions differs from that of the MSSM in two main aspects [23]: 



1. the LSP is no longer stable since it is not protected by symmetry, that is, it can 
decay within the detector; and 

2. it is possible to have single production of supersymmetric particles, since the final 
state is no longer restricted to be .R-parity even. 

We present the experimental status of SUSY particle search in imparity conserving case 
and i?-parity violating case separately. 



3.1 R-parity conserving case 
3.1.1 Sneutrinos and sleptons 

Once enough energy is available, the most efficient way to produce superpartners is to 
produce them in pair in high-energy leptonic or hadronic interactions. Sneutrinos, which 
are most likely invisible like neutrinos, can be produced in pair in the Z° boson decay at 



LEP. The partial decay width for sneutrinos is related to that for neutrinos by [24] 



T(Z°^v e V e ) = I/? 3 r(Z° ^i/^), £ = e,u or r, 

where ft = v~jc. The invisible decay width of the Z° measured at LEP is in good agreement 
with the standard model of N u =3 within errors. Therefore any additional contribution to 
the invisible decay width of the Z° is limited by 

invis. 
AT < 8.8MeV, 

at the 95% CL. This leads to the mass limit of m~ > 43 GeV if the three sneutrinos are 
degenerate in mass, while m~ > 40 GeV for one generation sneutrino p5| . 
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Events of the charged slepton pair production through the Z° decay 

Z° — > £ + £~ — > + 2 unobserved neutralinos 

are characterized by non-coplanar lepton pairs with large missing energy -momentum in 
the detector. The expected main backgrounds are Z° — > t^t (7),Z°^e + e £ + £ , Z° — > 
e + e~ + hadrons. The 95% CL lower mass limit on e, ju and t obtained at LEP is about 



45 GeV almost independent of slepton species 1 25] 



3.1.2 Gluinos and squarks 

Since gluinos and squarks are strongly interacting, they would be the SUSY particles with 
the largest cross sections at pp colliders. The branching ratios of gluinos and squarks 
decaying into various chargino and neutralino states depend on respective masses and 
mixing angles. In a scenario of a low mass gluino and massless photino (7 ) as the LSP, 
the decay modes become very simple. If m~ < m~ , the squark decays dominantly into gq 
and the main gluino decay is g — > qqj. If m~ < m~ , then the gluino decays dominantly 
into qq and the main squark decay is q — > qj. Since the photinos escape the detectors 
without any signature, the SUSY events would have two or more jets with a large amount 
of imbalanced transverse momenta. 

The limits on m~and m~, presented by the CDF experiment, are based on the compar- 
ison of the observed missing transverse energy (Ifix ) distribution with predictions for the 
standard model background and the QCD background plus SUSY contribution estimated 
by the ISAJET Monte Carlo samples p6| . Squark and gluino mass limits obtained at the 
CDF experiment , under the condition of SUSY parameters [/, = —250 GeV, tan/3 = 2 
and m#=500 GeV, are shown in Fig.l together with those of DO, LEP and UA1/UA2 
experiments (27j. A search for squarks and gluinos was made by DO in the three or more 
jets plus Ifix channel. The number of events observed in the data sample from an inte- 
grated luminosity of 13.4 ±1.6 pb" 1 was consistent with background. The limits on m~ 
and m~ plane are shown in Fig.l . For heavy squarks, a lower gluino mass limit of 146 
GeV was obtained, and for equal squark and gluino masses a mass limit of 205 GeV was 



obtained at the 95% CL [27]. 



3.1.3 Charginos and neutralinos 

Charginos are the mass eigenstates corresponding to the linear combination of winos and 
charged Higgsinos. Charginos have been searched for at LEP via their pair production 
and decays 

Z° -► WfWi, 
qqZ\,lv 'Zi, 

The mass limit at L3 comes from the line-shape measurement constraint : > 44 

GeV. This result is independent from the chargeno decays and the field contents(wino and 
charged Higgsino mixing) p8| . Searches for the signatures of acoplanar leptons, acoplanar 



e e 
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jets and isolated particles were combined at ALEPH. Irrespective of the field content and 
of the mass of the LSP, masses for chargeno below 45.2 GeV are excluded at the 95% 
CL|9|. 

Neutralinos, denoted by Z\, Z 2 , Z3, Z^ in order of increasing mass, are linear combina- 
tions of the photino, the zino, and the Higgsinos. The neutralino pair production at LEP 
is searched for via the following processes: 



€ re — > Z —> Z1Z1, ZiZ 2 , Z 2 Z 2 . 

In most cases, the main Z 2 decay mechanism is 

Z 2 -^ZiZ*^ Ztff or Z 2 -> Z17. 

Since the Z\ escapes undetected, the signature of the neutralino pair production would 
be excess in Y mv ^ s ' measurement. The signature of the above Z 2 decay processes are 
missing energy due to the undetected Z\ and one or two photons, two or four acollinear 
and acoplanar leptons, or one to four hadronic jets from the primary quarks [31|. 



From the absence of any candidate of Z\ and Z 2 at DELPHI [ 30 1 and L3[31|, limits 
on the probability that the Z° would decay into these channels were set. Combining 
these results with those from the Z° widths within the MSSM, one can caluculate, for 
every combination of M, fi and tan (3 values, the total contribution of Z° decays into 
neutralinos and charginos, and can deduce the excluded regions in the (M, fi) plane. The 
excluded regions for four values of tan (3, obtained by the L3 experiment through search 
for e + e~ — ► Z° — ► Z\Z 2 or Z 2 Z 2 are shown in Fig. 2 |ijf| . For moderate or high values 
of tan 0, a significant part of the accessible parameter space is excluded. All neutralino 
masses are functions of the parameters M, fj, and tan (3. Therefore, constraints on the 
MSSM parameter space translate into limits on these masses which are shown in Fig. 3 as 
a function of tan /3[]31|j, 

Data taken at the TOPAZ and VENUS with the TRISTAN e+e~ collider have been 
analyzed to study single photon events. The TOPAZ , based on data of an integrated 
luminosity of 164 pb _1 at the energy yfs = 59 GeV, observed 4 single photon candidates 
remained after event selection, which are consistent with the prediction of the standard 
model plus background |32]|. 



The VENUS , based on data of 225 pb _1 at \/s = 58 GeV, measured the single photon 
cross section , which is consistent with expectation. No anomalous signal has been observed 
leading to a lower limit on the mass of SUSY particles under assumption of radiative pair 
production of photinos. For massless photinos, the scalar electrons of the mass degenerate 



case are excluded below 51.9 GeV at the 90% CL [33]. Combining this value with other 



single photon experiments, the lower limit on the scalar electron was determined to be 
72.6 GeV at the 95% CL. Fig. 4 shows the excluded region on the (m~, m~) plane together 
with the results from other experiments compiled in Ref . [j33jl . 

Bounds on the gaugino parameters /i, tan (3 and M 2 have been calculated from exper- 
imental data[|J]. Fig. 5 shows the region in the (n,M 2 ) plane for tan/3 = 2 excluded by 
the experimental data on 

A lower bound on the mass of lighter chargino > 45 GeV, 
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B upper bound on the branching ratio of the visible neutralino mode Br(Z° — ► vis.) < 
5 x 10" 6 , 

C upper bound on the invisible width of the Z° , and 

D accepted gluino mass bound at CDF m~ >150 GeV (90% CL), 

where the hatched regions of each contour have been excluded. We see that the stringent 
bound comes from the constraint on the visible width of the Z° at LEP (B), or from the 
gluino search at CDF(D). 

3.1.4 Stop 

Among squarks, it may be that the stop is very light because of the high mass of the top. 
For a particular value of the mixing angle between tji and t^ , the lighter mass eigenstate 
t\ decouples from the Z° . Pair production of light stop has been extensively searched 
for by OPAL in e + e~ collisions at LEP ]35|]. The decay mode t\ — ► cZ\ was searched for 
in the data of the integrated luminosity of 69.1 pb -1 , which corresponds to 1.68 x 10 6 
produced Z° — ► qq events. With no t\ candidates, the region where more than 3.0 events 
are expected is excluded at 95% CL. Fig.6 shows the excluded region in the (9t,mj) plane 

for the case of the mass difference Am between t\ and Z\ larger than 5 GeV together 
with the region excluded by lower energy experiments [|35|] . The structure in the OPAL 
limits on 9t in Fig.6 is due to decoupling of the t\ from the Z® for 9t ~ 0.98. The OPAL 
excludes the existence of the t\ with a mass below 45.1 GeV at the 95% CL, where 9% of 
tji and ti is smaller than 0.85 or greater than 1.15, and Am is greater than 5 GeV. The 
exclusion regions in the (rar ,tti~ ) plane are shown in Fig. 7 for various 9t- The DELPHI 
has excluded significant regions in the (rar , m-x ) plane for different mixing angle 8t by 

assuming that t\ — > cZi, the most likely decay mode allowed by existing limits (3^]. The 
stop below 46 GeV has been excluded for any value of the mixing angle, except for a small 



window in 9t ~ 1 [ 30 [ . 



A search for the light stop has been carried out by the VENUS at TRISTAN [37]. A 
data sample of 210 pb -1 has been analyzed to find events with large acoplanar particle 
groups, assuming the two-body decay t\ — ► cZ\. The observed number of events was 
consistent with that expected from the known processes. They obtained the mass limits 
in the (mr , m~ ) plane at 95% CL, where t\ was excluded for the mass region from 
7.6 GeV to 28.0 GeV for massless LSP . The inclusive D*^ production cross section in 
two-photon processes was measured with the TOPAZ detector p8[| . The differential cross 
sections da{D* ± ) / dpT obtained were compared with theoretical predictions, such as those 
involving direct and resolved photon processes. A discussion was made on interpretation 
of the data of da^D^/dpr assuming stop pair production decaying into a charm quark 
and a neutralino Ml. 



A constraint from b — > 57 process to the MSSM has been derived in the light stop region 
It was pointed out that although some region in the parameter space is excluded from 
this process there remains a large parameter space where the amplitude of the b — > S7 
is suppressed due to partial cancellation between different diagrams. Stops as light as 20 
GeV are still viable from the b — > S7 constrainrt. It is also pointed out by Fukugita et 
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al. that the light stop with a mass as small as m^-~ 20 GeV is still an allowed possibility 
within the standard, minimal SUSY GUTs with the SUSY breaking induced by minimal 
supergravity [41]. 

Motivated by the fact that the stop may be considerably lighter than other squarks, 
H. Baer, J. Sender and X. Tata have reinvestigated its signals at the Fermilab Tevatron 
with simulations using ISAJET 7.07 under the assumption that it decays via t\ — > bW\ 
or t\ — ► cZ\. They have shown that experiments should be able to probe stop mass up to 
~ 100 GeV with an integrated luminosity of 100 pb _1 p^ ]. 

3.2 R-parity violating case 

3.2.1 Constraints from low-energy processes and neutrino physics 

We are particularly interested in the case in which the term X'^LiQjDf, in Eq. (4) in 
Sec. 2 is dominant, since this leads to resonant squark production through the s-channcl 
in the e-q subprocess at HERA |43|] . 

The term \' llk LiQiDi, (k = 1,2,3) contributes to the semileptonic decays of quarks, 
and the term X'^^xQ^D^^ (k = 1,2,3) contributes to forward-backward asymmetries 
measured in e + e _ collisions Existing limits for these coupling constants are given 
in [16]as a function of the mass of squarks. The term X' 1 j 1 LiQjDi, (j = 1,2,3) also 



contributes to atomic parity violation |16j]. The la bound for the coupling of this term is 

X' in < 0.26(m~j /lOOGeV), 

where the effect of the radiative corrections is taking into account^]. A general analysis 
of the constraints derived from neutrino physics on explicit i?-parity breaking in supersym- 
metric models is presented |44j|. The upper bounds for A^ and A^ - fe were obtained through 
(i)neutrino oscillations, (ii) neutrinoless double beta decay and (hi) neutrino cosmology 
and astrophysics, which were stringent in particular for lepton number violation involving 



the third generation [ 44 1. 



3.2.2 HERA 

ii-parity violating SUSY particles as well as leptoquarks, leptogluons and excited leptons 
are being searched for at the HERA ep collider by HI. In a data sample of ~ 0.5 pb _1 
no evidence was found for any squark production in the mass between 45 GeV up to 275 



GeV [45, 46]- In the super symmetric model, the single production of squark is allowed 
through imparity violating couplings. If the model is maximally imparity violating or if 
the photino mass is larger than the squark mass, these states are indistinguishable from 
leptoquarks. It implies the rejection limits for a search of squarks from ii-parity violating 
supersymmetry. The rejection limits on the coupling A' m are shown in Fig. 8 as a function 



of the mass of squark for various fixed photino masses [46]. 

The new particles are being searched for at HERA by ZEUS. No evidence for lepto- 
quark, leptogluon, squark, or excited electron production was found in a data sample of 
26 nb _1 |47j] . Limits on the production of squarks were determined for masses above 25 
GeV. Using the leptoquark limits, one can derive limits on the i?-parity violating coupling 
at the 98% CL as shown in Fig. 9. For electroweak coupling (^^ir/a), this yields lower 
limits on squark masses at the 95% CL : mj>168 GeV and mj~ > 92 GeV 
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4 Production processes at HERA 



4.1 Scalar electron production processes 

HERA is the world first ep collider, at which we can observe various inelastic processes, 
i.e., eq, eg, jg collisions and so on, in addition to the elastic processes ep — * Xp. For the 
purpose of the sparticle search, each collision processes could be useful for us to detect 
signatures from different kind of sparticles. In this subsection we discuss the scalar electron 
(selectron) production at HERA. 



4.1.1 R-parity conserving case 



The simplest sparticle production process at HERA would be the SUSY neutral current 
process (|| ||] 

ep -► eqX, (21) 

where e and q denote the selectron and squark, respectively. Even in such a simple 
process, e and q could be produces at the same time and we could get information of both 
masses of e and q from analyses of the process. This is a remarkable advantage of the 
ep collider HERA. Phenomenological analyses for the process have been given by many 
authors jl?], 5C] and we can find that the total cross section of the process will be cr~0.1 
150 GeV. Unfortunately, however, the TEVATRON experiment has 



pb for m~ + m~ ~ 
been excluded a light squark, m~~ 



150 GeV, by the negative search for the squark pair 
production pp — ► qq*X [26, 27], where they assume 5 flavor squarks are degenerated in 



mass. 



Therefore, we should say that the search for e and q by (21) will be not so hopefull. 
From the same reason the SUSY charged current process [48, O ] 

vqX (22) 



ep 



and the squark pair production via boson-gluon fusion [52| 

ep eqq* X (q / t) 



(23) 



will have too small cross section for us to extract from some large background processes. 
Note that the large lower mass bound from the TEVATRON could not be applicable to 
the stop. 



The SUSY bremsstrahlung process 



ep — > eZiqX 



(24) 



has also been analyzed by several authors [|53|, 54, 49] since this process has the dominant 
contribution from the t-channel photon exchange diagram which is not affected by the 
existence of heavy squark. However, the total cross section for this process will be rather 
small cj~0.01 pb even for the light selectron, m~~ 50 GeV, which has been almost excluded 
by the LEP. 
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Up to now we have discussed the inelastic processes only. On the other hand, the 
elastic selectron production 

ep — > eZ\p (25) 



has been known to have rather large and viable cross section [54, |55| . In particular, Lopez 
et al. |56| have pointed out recently that the Leading Proton Spectrometer at HERA 
would be efficient detector, which can measure momentum of the leading proton scattered 
in very forward region. 

The slepton - squark production in the NC and CC processes at LEP (g) LHC whose 



total available energies are larger than HERA is also discussed in Refs.[57, 58]. 



4.1.2 R-parity violating case 

If we consider the -R-parity breaking moel, the experimental signatures of sparticles will be 
completely different from those in the -R-parity conserving case as has been shown. First, 
the lightest sparticle will not be stable and will decay into the ordinary particles via the 
R-parity interaction. Second, some sparticles could be singly produced in eq sub-processes 
at HERA via the R-parity coupling A' in Eq.(||). 

In the former aspect, Dreiner et al. 5£] have shown that the SUSY neutral current 
process Eq.(|2~l|) will have remarkable signature and the mass reach at HERA will become 
about m~ + m~ ~ 200 GeV even for very small R-breaking couplings A ~ 10 . 

The extensive study for the latter situation has been given by Butterworth et al. fi~4| 



and by ourself [13, 60, 61]. The single production of the SUSY partnars of light quarks 



has been discussed in Refs.|14j, |15[| . On the other hand, we have published some papers 
in which considered the single stop production at HERA. Our works will be summarized 
in next subsection. Butterworth et al. have analyzed extensively the single production of 
squarks with RB couplings in the first and second generation H]. For example, 

e~u -> d (26) 

could be possible and the d-squark will decay via .R-conserving couplings 

d^dZ x . (27) 

Moreover, the lightest neutralino Z\ will decay into ude~ or ude + through the .R-breaking 
interaction proportional to A'. Here the interesting point is that the decay product could 
contain not only the electron but also the positron. This is because the neutralino is the 
Majorana fermion. In fact, the positron signature from the electron-proton collision will 
be distinguishable from almost standard background. In their model of cascade decays of 
the squarks to the LSP, the squark mass reach was formed to be at HERA 

m~ ~ 270GeV for A' ~ 0.08 

and the reach in the Yukawa coupling is 

A' ~ 5.3xl0~ 3 for m~~ lOOGeV. 
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4.2 Stop production processes 
4.2.1 Boson-gluon fusion 



The purpose of the present subsection is to discuss the production of the light stop at 
HERA in the framework of the MSSM with the i?-parity conservation. The possible exis- 
tence of the light stop t\ has been discussed by some people [H], ^2], 34, |4l| . The most 
promissing production process in ep collisions will be the boson-gluon fusion(BGF)|33|, |5^, 
6§ : 

e ± p -> e ± t 1 t* 1 X. (28) 

Its Feynman diagrams are depicted in Fig. 10. Although it will be shown that the Weizsacker 
- Williams approximation (WWA) is very appropriate to our purpose, we first carry out 
the exact calculations 



da = ^G{r ) ,s) J di 1 dPS^x^Y.\ M a + M b + M c \ 2 . 



(29) 



spin 



Here r\ denotes the momentum fraction of gluon in the proton and s = (pj + Pf) 2 - The 
gluon distribution function G(rj,s) is given by the set 1 of Ref.[64| and the three-body 
phase space volume is expressed as 



dPS {3) 



1 



{ . )%3 -e( S + 2m 2 e ~ ^m 2 e (s + m 2 ) - W 2 )dydQ 2 dz^- 



with 



W 



l = + m vf 



m r 



Z =p- p f /p ■ q, 

cos< ^(PXle)-(pX P/ ) 

|pxl e | • |pxp/| 

The invariant amplitudes corresponding to Fig. 10 are, respectively, given by 

ie 2 g s (q - 2p f ,) fl (p - 2p f ) L 



M a 

Mr ■■ 



2 

m£ 



t 



Q 2 



ie 2 9s (<? ~ 2 Pf)v(P ~ 2 P'f)v 
2ie 2 g t 



2 

m£ 
h 



Q 2 



with i = (pj — p) 2 ,u = {pf — p) 2 and 



rp/lU 



+ o 



Q 2 



h Q 2 + m| 



(AlP L + A%P R )Ul e ) e u (p) 



G 



cos 2 9 t 



| sin 2 9w 



2 cos 6 W sin 6w 



i e L = -(cotOw — tan^V), 



Pl 

R 



1 



tan 6^, 
lT7s)- 



(30) 



(31) 
(32) 
(33) 
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It is worth mentioning that the constant C~ representing the strength of t\t\Z cou- 
pling depends sensitively on the mixing angle 0% [22]. In particular, Or = if 9t = 
cos _1 (-^ sin 9\y) ~ 0.98. On the other hand, i*ii7 and t\t\g couplings do not depend on 
Of Therefore, the Z boson contribution to the cross section depends on 9t, while that of 
the photon is independent from 9%. 

It is interesting to compare our exact tree level calculation with the Weizsacker- 
Williams approximation (WWA) ; 

^=J dyP(y) J dr,G(r,,s')^, (34) 

where s' = yr]s. The WWA factorizes the cross section of the process, shown in Fig. 10, 
into the probability for emitting photon from lepton ; 



' mill 



and the BGF cross section involving real photon ; 
da 4 1 



dz 9 naas s' 3 z 2 (l - z) 2 L """*i 

The numerical calculations have been performed by using the program packages BASES 
[65]. The standard model parameters are taken as a = 1/137, sin 2 9\y = 0.23 and mz = 



[2mA _ 2mf s'z(l - z) + s' 2 z 2 (l - z) 2 }. (36) 



91.1 GeV. In Fig. 11 we show the stop mass dependence of the total cross sections a in 
the exact tree level calculation and in the WWA. It is found that the 7 contribution is 
much larger than the Z boson contribution, therefore the WWA is a good approximation. 
As mentioned already the photon contribution does not depend on the stop mixing angle 
9f. The total cross section, consequently, is insensitive to 9t- Throughout the calculations 
in Fig. 11, we have adopted the lower cut for Q 2 as Q 2 > 5 GeV 2 in order to make 
detectable the scattered electron. Because of the photon dominance, the total cross section 
significantly increases when smaller Q 2 cuts are adopted. If Q 2 is cut at the kinematical 
limit ; Q 2 > m 2 -^—^, the total cross section is about four times as large as that with cut 

Q 2 > 5 GeV 2 . If we adopt Q 2 cut Q 2 > m 2 -^r , the detectable cross section turns out to 
be a > 0.1pb in the case of mr < 50 GeV. 



4.2.2 Resonance production 

In the present subsection, we shall concentrate our discussions to -R-parity breaking (RB) 
process. We start from a coupling of the stop t\ 

C-int = A' 131 cos 9 t (tidP L e + t\eP R d) (37) 

originated from the RB superpotential Eq.(0). Here Pl r denore left and right handed 
chiral projection operators. The coupling Eq. ( |20| ) is the most suitable for the ep collider 
experiments at HERA, because the stop will be produced through the s-channel in the 
e-d subprocess as shown in Fig. 12 for example, 

e T p -» (hX) -» e T qX. (38) 
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For simplicity, we will assume A' 131 to be only non-zero coupling parameter in the following. 
The upper bound on the strength of coupling has been investigated through the low-energy 
experiments ]l6| and the neutrino physics §01]. The most stringent bound A' 131 ~ 0.3 comes 
from the atomic parity violation experiment [ill ]. 

As mentioned in 2.3.3, only two-body decay modes (a), (b) and (h) are targets of our 
study. The formulae of the decay width for each mode are respectively given by 



tZi) 



a 

ti 



A 2 (m~ , to? , to- ) 



(\F L \ 2 + \F R \ 2 )(r4-ml 



4m t m~Re(F R F* L ) 



Fr 



m t NH cos 6 t 
2m\v sin 9y/ sin f3 

1/2 - e u sin 2 9 W 



+ e u (N% - tan W N'* A ) sin Q u 

m t iV| 4 sin 6 t 



cos 9w sin 9\\r 



N' i2 ) cos i 



2niy/ sin 8\y sin j3 ' 



(39) 
(40) 

(41) 



r(ti - 6W fc ) = — (mf , to 2 , to 2 - ) 

|G L | 2 + \G R \ 2 )(ml - m 2 - m 2 ~) - Am b m~Re(G R Gl)] , (42) 
_ m b U£ 2 cos6 t 

l = — -m ( 43 ) 

V zmvy cos p 

^ -t/- /} i ^Vfeasin^ 

= Vfci cos 6> t + -= — , (44) 

V 2mw sin p 



and 



r(ti ->ed) = ^cos 2 ^, (45) 

where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. NL, Vki and Uki respectively stand for 
the neutralino and chargino mixing angles || . The mixing angles as well as masses of the 
neutralinos m~ and the charginos m~ are determined from the basic parameters in the 
MSSM(//, tan/3, M 2 ). As seen from Fig. 13 the branching ratio of the stop depends on the 
stop mass. If the stop is heavy enough, i.e., rrij^ > mt + m~ or > to& + and the 

RB coupling is comparable with the gauge or Yukawa coupling A / 1 2 1 /4-7r ~ a, at there is 
a parameter region where BR(ti — ► tZj) or BR(t\ — ► bW^) competes with BR(t\ — > ed). 
For light t\, however, the decay t\ — > ed predominates over other decay channels and we 
have BR(h -» ed) ~ 100%. 



(A) The case of light stop : mr < rat + m~ or ra~ < mi, + TO 



Wk 



Assuming the stop has BR(t\ — > ed) ~ 100%, we can calculate the inclusive differential 
cross section for e^p — > e^gA with polarized e 1 * 1 beams(see Fig. 12) as follows: 
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^2 Kr] = Q 2 ) E Ti{e% R q) + q(x, Q 2 ) £ Jtfe^g)], (46) 

where the coefficients Ti(e±L,Ftq) are represented in the Appendix. For e^d — > ti followed 
by the decay ti — > e^d we have 



where the decay width T(t\ — ► ed) is given by (45). From Eq.(47) we can expect a clean 
signal of the stop as a sharp resonance peak in the distribution of the Bjorken parameter 
x. The position peak corresponds to x = rr& /s for fixed Q 2 . Figure 14 clearly shows 
the resonance behavior of the production process. The background curve represents the 
predicton of SM. It is obviously seen from Fig. 14 that the lower Q 2 cuts would be very 
efficient to suppress the background expected from SM, since the s-channel resonance 
contribution is independent of Q 2 . It would be worthy of remarking that the similar peak 



could be expected in the leptoquark production at HERA[p6|. We should point out that 
the stop with the RB couplings will be discriminated from most of the leptoquarks by its 
distinctive properties: (1) the x peak originated from the stop would exist only in the NC 
(not exist in the CC) process due to no RB stop couplings to neutrinos, (2) e + beams are 
more favorable than e~ beams as will be mentioned later. One of the leptoquarks 51/2 
with the charge Q = — | will give the same signature as the RB stop, if the stop decays 
into the electron and ci-quark with BR(t\ — » ed) ~ 100%. Q The event rate depends on 
not only the RB coupling strength A' 131 but also the kind of beams. Figure 15 shows the 
y distribution at fixed x for e~ and e + beams. We can see that the e + beam is more 
efficient than the e~ one to distinguish the stop signal from the SM background. This 
can be understood from the fact that the e + collides with valence (i-quark in the proton, 
while the e~ does only with sea d-quarks. The difference of their distribution functions 
is naturally reflected in the cross sections in Fig. 15. It is hopefully expected that the 
longitudinally polarized e + and e~ beams will soon be available at HERA. They could 
be advantageous to suppress the SM background. In Fig. 16 we show the y{= Q 2 / s) 
dependence of the asymmetries defined by 



Cr = da(e^)/dxdy - da(e R )/dxdy 
da(e R )/dxdy + da(e R )/dxdy 



and 



A _ da{e L )/dxdy - da(e R )/dxdy ^ 
da(e2)/dxdy + da(e R )/dxdy 

It will be seen that the longitudinally polarized beams will be useful to identify the stop 
signal for a reasonable range of y in Cr . Finaly we show the searchable parameter region 

1 H1 collabaration at HERA has given the lower mass bound m~ ~ 98 GeV on the RB stop from the 
negative result for the leptoquark S1/2 search at 95% CL or A' 131 = 0.3 E3]. 
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in (A' 131 , m,~ ) plane at HERA. In Fig. 17 the shaded region is experimentally excluded from 
the atomic parity violation experiment [16]. The area inside the solid contour is accessible 
at HERA whose production rate is more than ten signal events above the SM background 
with Q 2 > 10 3 GeV 2 in 100 pb _1 running. From Fig. 17 it is seen that the stop mass 
reach is 

to- ~ 200(270)GeV 
for e~(e + ) beams at the RB coupling 

A'isi - 0.1 



(B) The case of heavy stop: > mt + m~ or m~ > m^ + 

If the mass of the stop is heavy enough, various decay channels as (a) to(g) compete 
with [h) and a sharp peak at x = mj 2 /2 can no longer be expected. For the case of 

BR(t\ —* ed) <C 100%, we should take into account of the processes 

ep -> tZiX (50) 

and 

ep -» bW k X (51) 

as mentioned before. Their Feynman diagrams are shown in Fig. 18. Here we should 
consider the virtual contributions of the selectron, sneutrino and d-squark with the same 
RB couplings constantsA^ 31 . The differential cross sections are given by 



da , - . aX 



a 

131 



(u — m 2 )(u — to-) (i — m 2 )(i — m~) 



(ep - tZ t X) = | F ~| 2 _„ - + \F4 



dxdQ 2 ^ ■ 8§ 2 L'*ei (u-m 2 ) 2 1 dl (i-m 2 -) 2 

V e L J ^ d R J 



(s — mi ) 2 — mi Ti 



+2Re(F 7 FZ 



tu — mi mi- 



(u — TO- ) (t — TO- ) 
e L d R 



2 cos 2 9 t s(s - mi ) . 

-Re[F^{F R u + F L m t m^ 



((s — mi) 2 + miTi)(u — mi- ) 
2 cos 2 9ts(s — mi ) 

Re{FUF R t + F L m t m~)\ , (52) 



((§ — mi ) 2 + mi- r~ )(u — to- ) 
where s = xs, t = —Q 2 and 



77 - at' 1/2 + e e sin 2 9w AT , 

F~= e e N a — N i2 , 

cos (?w sm 0w 

F~= eeN^-edtanOwN^. 
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da 



dxdQ' 



: (ep^bW k X) 



16s 2 sin 9w 



\Vi: 



(u — m?)(u — m 2 ~ ) 
(u — mi) 2 



+- 



cos o t s 



(s — mi ) 2 — mi T~ 
v W ti ti 



{\G L \ 2 + \G R \ 2 ){s-ml 



m- 



4m b m^ Re(G R G* L 



2 cos 2 O t s(s — m~ ) 



(s — mi ) 2 + mi Ti 



[u — mA) 



Re{V*i{GjiU + G^m^m 



(53) 



Also from Fig. 19 we can infer that e + beams are more efficient than e ones. It is 
expected that the detectable cross sections a ~ 0.1 pb for heavy stop with mass m~ ~ 250 
GeV for e + beams. As far as e~ beams are concerned only e~p — > bW^X would be 
detectable for mr ~ 170 GeV. In our model the LSP, the lightest neutralino Z\ will decay 
into i?-even particles via only non-zero RB coupling A 131 . A typical decay chain will be 
(see Fig. 20) 

ep -> tZiX. -» (bW){bdu)X -» (b(iv))(bdv)X. (54) 

The chargino W\ will also decay into i?-even particles through Z\ like 

ep -» bW\X -» (MvZ^X -» b(£u(bdu))X. (55) 

In both processes ( ^0|) and (|5l|) 2b-jets+jet+lepton+missing transverse momentum 
would emerge as final products to be detectable. Figure 21 shows the Monte Carlo events 
for the transverse momentum distribution of scattered muon from the process (]55| ) under 
the condition of the integrated luminosity C = 300 pb -1 for e~p collisions at HERA. The 



branching ratio BR(W\ — ► v\iZ\) is assumed to be ^ [67] . Shown in Fig. 21 are distributions 
for m~ = 150 GeV and rrv^ =100 GeV together with background(short-dashed line) of 
muon events coming from charged current process(CC) e~p — > uqX and VF-gluon fusion 



process (WGF) eTp -> vscX, vbcX. The generators LEPTO )6J] and AROMA [|69| with 
JETSET [frcfl have , respectively, been used for CC and WGF. From Fig. 19 we can see 
that ep — * bW\X will have a reasonably feasible cross section to which the stop contributes 
from the s-channel for m~ ~100 GeV. This is also the case for ep — ► tZ\X since the same 

final states are realized. Thus, in both processes, ep — > tZiX&nd ep — * bW\X a possible 
typical signature of the stop production would be 2b-jets+jet+lepton+missing transverse 
momentum ^j* owing to the LSP decay via RB coupling. One of the signals to be detected 
at HERA is characterized by the high Pt spectrum of muons where the lower Pt cut makes 
the event distinctive from its background. Although it is difficult to discriminate the stop 
from one of the leptoquark Si/ 2 for light enough stop m~ < mt + m~ ,mj^ < m^ + m ^ k -, 
a characteristic signature of the heavy stop with mass m~ > m t + m~ , > + m^ 

could clearly be distinguished from the leptoquark Si/ 2 - 

Recent observation of a single muon of high Pt event in e + p — ► fi + X by HI jTl] ] gives 
us a favourable evidence for our arguments on the stop. As seen from Fig. 19 a few events 
would be expected for m~ = 150 GeV and present luminosity L ~ 3 pb -1 . 
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5 Concluding remarks 



We have discussed SUSY particle searches at HERA by paying our particular attention 
upon the single production of stop in the framework of the MSSM with i?-parity breaking 
interactions. Observation of the moderately light stop predicted by the model would be a 
nice target of HERA experiments so far the RB coupling constant A' 131 ~ 0.3. 

In the case of BR(ti — > ed) ~ 100% corresponding to < m t + or m~ < 
nib + m y^ 1 > the stop produced via neutral currents can clearly be seen as a resonance peak 
at the Bjorken parameter x distribution. Heavier stop, m~ > nit+m-x or m~ > m^+m- , 
are produced through charged currents emerge from the following chain of decay processes 

ep -► tZtX -► {bW){bdu)X -► {b{iv))(bdv)X. (56) 

The chargino Wi will also decay into i?-even particles through Z\ like 

ep -» bW\X (MvZ{)X b{lv{bdv))X. (57) 

The evidence for the existence of the stop would emerge from the careful identification of 
26-jets + jet + £ + ^V- Possible backgrounds come from e~p — > and — > z/scA, 
vbcX. The transverse momentum distribution of scattering muon will distinguish the stop 
events more clearly from backgrounds. Extensive Monte Carlo events with a variety of 
SUSY parameter sets would highly be desirable. Since the production cross section of the 
stop is unfortunately rather low the high luminosity is very much favourable. 

Until LEPII or Next Liner Collider will be built we are sure that HERA would play 
a unique role to open the "SUSY world" through our stop with mass of 100 to 300 GeV. 
TEVATRON groups are also enthusiastic about sparticle search at high energies. We are 
very much expecting HERA and TEVATRON complementarily step forward to enter into 
the novel region of particle physics by making full use of their specialities of detectors. 
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Appendix 



The analytic expression for the cross section Eq.(filf) is given as follows. 
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(60) 
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(62) 

(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 
(71) 

(72) 
(73) 



A\A R \ (74) 

(75) 
(76) 

(77) 



22 



for e beams. The formula for e + beams can be obtained by the following replacement in 
the above formula for the e~ beams ; 

e L,R e R,L ( 78 ) 

q - q (79) 

q - g (80) 
Here, denote the electromagnetic charge of matter fermion /, and 

A f = Tj-efsm 2 e w ^ 

L cos 9\y sin 0^ 

A f R = e f tan6 w , (82) 

where are the third component of isospin and 0w is the Weinberg angle. 
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Figure Captions 



Figure 1: DO, LEP and UA1/UA2 squark and gluino mass limits as a function of squark 
and gluino massp 7 / 



Figure 2: The excluded regions of the MSSM parameter space at 95% CL as a function 
of the parameters M and /i. The kinematical limit corresponds to the sum of the 
lightest and next to lightest neutralino masses being equal to the center-of-mass 
energy. The exclusion coming from the lineshape measurement goes beyond the 
kinematic limit of the direct search |M| . 

Figure 3: Neutralino lower mass limits(95% CL). The lines correspond to the lower mass 
limit for the different neutralinos, where Xi x'j x" an d x'" i n this figure represent 
Z\, Z%, Z% and Z4 in the text, respectively. The regions below the lines are excluded 



m 



Figure 4: Excluded region on the (m~, m~) plane. VENUS (90% CL, solid line), ASP 
(90% CL, dashed line), CELLO (90% CL, dotted line), ALEPH (95% CL, dot-dashed 
line) and combined single photon result (90% CL, solid line) 

Figure 5: The region in the (//, M2 ) plane for tan (3 = 2 excluded by the experimental 
data on 

A lower bound on the mass of lighter chargino m^. > 45 GeV, 

B upper bound on the branching ratio of the visible neutralino mode Br(Z° — ► 
visi.) < 5 x 10 -6 , 

C upper bound on the invisible width of the Z° , and 

D accepted gluino mass bound at CDF m~ > 150 GeV (90% CL), 

where the hatched regions of each contour have been excluded p4[. 



Figure 6: The excluded region in the (Oti'&rnix) , fn^ ) plane at 95% CL, for the case of 
the mass difference Am between t\ and Z\ larger than 5 GeV. The region excluded 
from the limit on the Z° decay width (Ar < 26 GeV at 95% CL) and limits from 
previous publications are also shown p5|. 



Figure 7: The excluded regions in the (m~ ,m~ ) plane at 95% CL, where the mixing 
angle is assumed to be Ot^&mix) < 0.85 or > 1.15 rad (shaded area), and Oti'&mix) < 
0.97 or > 0.99 rad (hatched area). The dashed curve shows the contour of the limit 
from the previonus publication |]35||. 



Figure 8: Rejection limits at the 95% CL for the coupling A' m as a function of the squark 
mass for various fixed photino masses. Regions above the curves are excluded. The 



limits combine all charged and neutral decays of the d and u [46]. 



Figure 9: The 95% CL upper limits on the couplings of squarks versus mass for (a) the d 
(in order of decreasing coupling limits based on the NC, CC and combined samples) 
and (b) the u and i using the NC samples only. The sensitivity of the result to the 



photino mass is also shown [47]. 
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Figure 10: Feynman diagrams for e g — ► e [52|. 

Figure 11: Total cross-sections for e~p — > e~iiifX with cut Q 2 > 5 GeV 2 . Solid line, 
dashed line and dotted line correspond to exact calculation including both 7 and Z, 
including only 7 and WWA, respectively [62]. 

Figure 12: Feynman diagrams for sub-processes e^q — ► e^q 

Figure 13: m~ dependence of branching ratio of stop. We take mt= 135 GeV, tan/3=2, 
9 t =1.0, \'i 31 =0.1 and (M 2 (GeV), /z (GeV)) = (50, -100) for (a) and (100, -50) for 
(b) 0. 

Figure 14: x distribution at fixed Q 2 . Adopted parameters are m-r =200 GeV, A' 131 =0.25 
and 6 t =0.0 @. 

Figure 15: y distributions at fixed x using the electron and the positron, x is fixed at 
0.2. Adopted parameters are m~ =140 GeV and 9t=ir/4 p0| . 

Figure 16: y distribution at fixed x of differential asymmetries Cr (a) and A e - (b). x is 
fixed at 0.2. Adopted parameters are m~ =140 GeV and 6^=0 [60| ]. 

Figure 17: Searchable parameter region at HERA in (X' 131 ,mt). Kinematical cut is Q 2 
> 10 3 GeV 2 IBCH. 



Figure 18: Feynman diagrams for sub-processes ep — ► tZ{ and ep — > |61]. 

Figure 19: Stop mass dependence of total cross section. We take m; =135 GeV, 6% = 
1.0, tan/? = 2, A' 131 = 0.1, M 2 = 100 GeV, m~ t = 200 GeV, m g - = 300 GeV and ft = 
—50 GeV. Solid, short-dashed, dotted and dashed lines correspond to e~p — > bWiX, 
e + p — > bW^~X, e~p — > tZ\X and e + p — > tZ\X , respectively |]6l[ . 

Figure 20: Feynman diagrams for LSP decay |61| . 

Figure 21: Monte Carlo events for transverse momentum distribution of scattered muon 
from e~p — > bWiX (solid lines) together with backgrounds CC and WGF processes 
(short-dashed line). We take m £ -=200 GeV, m g -=300 GeV , mt=135 GeV, t =l.O, 
tan/3=2, A' 131 =0.1, M 2 =100 GeV, ^=-50 GeV and integrated luminosity L=300 
pb- 1 ||. 
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Table I Particle content 



Superfields 


Spin-0 


Spin- 1/2 


Spin-1 


Color 


T 


Y 


B 


L 


Q 




Ql 




3 


1 
2 


+1 


+1 





u c 


UR 


K)l 




3* 





2 

3 


1 

3 





D c 


d-R 


(d c ) L 




3* 





+\ 


1 

3 





L 


£l 


i L 




1 


1 
2 


1 

2 





+1 


E c 


&R 


(e c ) L 




1 





+1 





-1 


G 




9 


9 


8 














W 




W 


W 


1 


1 











B 




B 


B 


1 














m 


Hi 


H\L 




1 


1 

2 


1 

2 








H 2 


H 2 


H 2 l 




1 


1 

2 


+\ 









30 



